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ABSTRACT 

We searched for evidence of the reddening of background QSO spectra from the Sloan 
Digital Sky Survey (SDSS) due to dust in intervening damped Lyman a (DLA) sys- 
tems. We utilise the two Data Releases DR5 and DR7 to arrive at sample sizes of 475 
(DR5) and 676 (DR7) absorbers, based on two different published lists of SDSS DLAs. 
Both samples span roughly the redshift range of 2.2 < Zabs < 5.2, with a mean of z^3.0, 
and the majority of the DLAs (75%) below z=3.3. We construct geometric mean spec- 
tra in the absorber restframes ranging from 1240 to ~ 2800 A, and composite spectra 
of samples matching the 'DLA' QSOs in i band magnitude and emission redshift Ze^, 
but without absorption lines. By comparing the slopes of these composite spectra with 
their matched counterparts, we find no sign of reddening in the ensemble of the ab- 
sorbers from these samples. Owing to both the unprecedently large sizes of the DLA 
samples themselves and the non-DLA SDSS QSO sample, from which we can draw our 
matching spectra, we can place very tight limits for this non-detection (<E(B-V)> 
=-0.0013± 0.0025 (DR5) and <E(B-V)> =-0.0017±0.0022 (DR7). Interestingly, when 
applying our technique to the samples of York et. al. (2006), vandenBerk et al. (2008) 
(intervening and intrinsic Mgll absorbers) and the smaller DLA-subsample and pool 
of comparison QSOs of Vladilo et al. (2008), we do recover their results, i.e. detect the 
same amount of reddening as these authors do. Furthermore, we have tested whether 
subsamples of our large sample in categories involving the absorbers (HI column den- 
sities, presence or absence of accompanying metal absorption, absorber redshift) or 
the background quasars (emission redshift, brightness) do reveal dust extinction, but 
found no trends. These results are at odds with both detections of dust reddening 
from previous studies, and also with expectations from observations of high-redshift 
galaxies. 
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1 INTRODUCTION 

Detecting galaxies at high-redshift provides a powerful 
observational way to probe structure formation in the 
early Universe. There have now been several hundreds 
of galaxies f o und using the Lyman-break technique (e.g. 
ISteidel et al.l l|l999l ). Although very successful, this method 
is limited to luminous objects. A complementary technique 
to select high-redshift galaxies is the detection of systems in 
absorption in the spectra of background quasars. Since their 
selection is independent of luminosity and morphology, 
the study of galaxy evolution through these systems cir- 
cumvents the biases inherent to traditional emission studies. 
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Damped Lyman-a systems (hereafter DLAs) have hy- 
drogen column densities, log N(HI) > 20.3 and are thought 
to be associated with parts o f high-re dshift galaxi es, just 
as th ey are at lowest redshift (|Rao et al. 2003; Che n et al.l 
l2005h . From our knowledge of the local Universe we expect 
that high-redshift galaxies, especially those forming stars, 
should contain some amount of dust. Therefore, the pos- 
sibility that DLAs are dusty objects is currently hotly de- 
bated. Indeed, smoothe d-particle hydrodynamic simulations 
(|Nagamine et al.ll2004l 'l predict that DLAs should be 1/3'^'* 
solar at z = 2.5 and even more metal-rich toward lower red- 
shifts whilst current observations measure typically a metal- 
licity around 1/10* solar. Dust is therefore invoked to recon- 
cile results from the simulations with current observations. 
In addition, one can estimate the dust content of individual 
DLA from the ratio of depleted onto undepleted elements it 
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contains. It is assumed that a fraction of the elements is not 
dete cted in the gas phase because it is locked into dust grains 
(e.g. IPettini et al.l l|2000h . Such studies have shown that the 
dust depletion levels observed in DLAs are small, lower than 
those measured in the Galactic interstellar-medium clouds 
with similar N(HI) column densities. In fact, the observa- 
tions resemble_th«dustde^^ halo warm 
gas (|Kulkarni et al.lll997l : ISavaglio et al.ll2000l ). 

Another more direct way to probe the dust content of 
DLAs is to look for a 2175 A feature (the so-called "dust 
bump") in the quasar spectrum. This feature is known to 
be more predominant in the Milky Way than in the Small 



detections (|Juralll977 
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l|2004 ) claimed the detections of bumps in three ind i vidua l 
Mgll absorbers at z=1.5, while lJunkkarinen et al] (|2004l ) 
have detected a clear bump in a DLA at z=0.524 toward a 
BL Lac at z=0.94. Other searches remained unsuccessful, 
indicating that most DLAs do not show a predominant 2175 
A feature. A promising tool for finding the signatures of dust 
are searc hes for the strong silica te absorption features in the 
infrared. iKulkarni et"all |2007l '). using Spitzer data, repor ; 
on the first indication of 9.7 pm silicate absorption in a 
DLA at 7.abs ~ 0.52, and find it to be stronger than expected 
from extrapolation of the relation from Galactic diffuse ISM. 

A consequence of DLAs being dusty is that it might 
prevent the inclusion of their background quasar in colour- 
colour selected samples. In order to overcome such biases, 
radio- selected quas ar samples have been built recently: 
while I2OOII l|EllisorJ ) did not discover a missing population, 
a more recent survey finds 4 fields (out of ~60) where the 
radio-flux selected quasa rs have no optical c ounterparts 
down to very faint limits jjorgenson et al.ll2006h . These ob- 
jects are the potential candidates of dusty DLAs obscuring 
background quasars. 

In this paper, we investigate yet another possible signa- 
ture of the dust content of DLAs. Indeed, dust in a quasar 
absorber will re dden the colour of its background quasar. 
iPei et all (|l99lD were the first to look for reddening of the 
quasar continua in a sample containing DLAs and detect 
reddening at 99.999% level. More recently. iMurphv fc Lisk3 
l|2004l) have used 79 DLAs from SDSS Data Release 2 to 
compare the slope of power laws fitting quasar continua with 
and without DLAs. These authors found a colour excess due 
to Small Magellanic Cloud type of dust of reddening <E(B- 
V)> <0.01 mag. Using colour excess measurements for SDSS 
sightli nes containing DLAs and matching sightlines without 
them, IVladilo et al.l l|2008l ) report a detection of reddening 
in a sample of DR5 SDSS QSOs towards sightlines contain- 
ing DLAs at the mean level of <E(B-V)> ~6.3xl0"^. They 
derive their estimate by comparing colours of QSOs with in- 
tervening DLAs to the mean colour of a set of QSO of similar 
emission redshift and brightness, and bootstrapping the re- 
sulting broad distribution of colour excesses to ascertain the 
slightly redder colours of the DLA QSOs. A byproduct of the 
iPontzen fc Pettinil (|2009l ) work, who performed a Bayesian 
analysis on a variety of optical and radio datasets to deter- 
mine the loss fraction of QSOs due to dust obscuration, is an 



estimate of the probability density distribution for dust red- 
dening in optical spectra, that also peaks around the same 
value for <E(B-V)> . Here, we build composite spectra of 
quasar spectra with DLAs from SDSS Data Releases 5 and 
7, and compare them with reference composite spectra built 
to match the original spectra in magnitudes and redshifts. 

The organisation of the paper is as follows : after the 
introduction, we describe the selection of the DLA samples 
in section 2, followed by the construction of the matching 
non-absorber samples, and the procedure for obtaining 
composite spectra in section 3, where we also presents the 
method for deriving extinction as well as the results for the 
full samples. After discussing subsample analyses (section 
4), we summarise and conclude in section 5. 



2 DLA SAMPLE SELECTION 

We start the selection of our samples from two lists of SDSS 

sur veys for DLAs. 

iProchaska fc Wolfd l|2009l) searched in Data Release 
(DR) 5, following a parti ally automated procedure tested 
on earlier data relea ses (|Prochaska fc Herbert-Fortl |2004| : 
IProchaska et al.ll2005h . Their cleanest ('statistical') sample 
contains 738 DLAs. We have chosen not to include objects 
that are listed in their non-statistical sample, as these 
can be towards QSOs with strong intrinsic absorption, or 
have Zabs ~ zqso- Broad absorption line QSO (BALs) are 
known to be more strongly reddened in the rest-frame UV 
than non-BAL QSOs, and hence this list of DLAs is then 
cross-correlated with the BAL catalogue of iGibson et al.l 
(20091), in order to avoid such BALs entering our analysis. 
Furthermore, we have r estricted ourselves to t hose s ightlines 
towards QSOs in the IProchaska fc Wolfd (|2009l ) sample 
that only exhibit one single DLA. Thus, we exclude the 59 
respectively 5 cases, where two or three DLAs are found 
along the same line of sight. This leaves us with a sample 
of 526 SDSS DR5 QSO sighthnes with securely identified 
DLA well suited for the analysis of their dust content. 



iNoterdaeme et al] (|2009l ) have surveyed DR7, and 
present a list of 1426 strong HI absorbers at redshifts 2.15^ 
z ^ 5.2, of which 937 systems have log N(HI ^ 20.3. Apply- 
ing the same criteria for filtering out potentially problematic 
sightlines as above, we arrive at 731 SDSS DR7 QSO sight- 
lines for this DLA sampleQ 

While both lists have a substantial overlap (371 of 
Prochaska's sources appear also in the Noterdaeme list, al- 
beit 31 of these are estimated there to fall below the log N(HI 
^ 20.3 criterion), the two methods for finding DLAs differ 
significantly enough to warrant keeping both datasets sepa- 
rate. It should, however, be kept in mind that the two sam- 
ples are not independent from each other. Note also that the 
SDSS treatment of spectroscopic calibration has improved 



^ The B AL catalogue can only ex clude objects that were already 
in DR5. iNoterdaeme et al.l ll2009h make an effort to clean their 
samples of such sigh tlines, but we nonet heless found 44 entries 
with a match in the I Gibson et al. I ll2009h catalogue. Hence, we 
also visually inspected our final sample of new DR7 sightlines to 
exclude obvious BALs, but found none. 



Dust Reddening in SDSS Spectra with Intervening DLA 3 



from DR5 to DR7, especially the spectrophotometric algo- 
rithms. Therefore, we have decided to strictly use only the 
versions of the spetra that have been processed by the DR7 
pipeline. 

Figure [T] shows the distribution of both samples in QSO 
emission redshift and the i band magnitudes. In addition, 
we show the underlying SDSS QSO lists (cleaned of BALs 
and the DLAs from both lists, as detailed below) from 
which we can draw our matching spectra. The insert shows 
the histogram of the absorber redshifts for the two samples. 
The vast majority (> 75%) of the DLAs are in the redshift 
range ~2.2 ^ Zats < 3.3. 



3 PROCEDURE FOR DETERMINING 
EXTINCTION 

There are a variety of different methods to de- 
rive estimates for the reddening of each QSO spec- 
trum , or th e average reddening of the whole sample 
fe.g. IPei et al. (1991); Ellison ct al. (20 051'): iRichards et all 
(|2003l ): iMurphv fc Liskel (|2004l ): IVladilo et al.l (|2008| ). We 
follow closely the sugge stion of lYork et al.l ()2006h and 
IVanden Berk et aP lj2008l ) , who have analysed the dust red- 
dening in Mgll absorbers towards SDSS QSOs. The basic 
principle is to define a sample of QSOs, where for each DLA 
quasar in question a close match in both emission redshift 
and i-band magnitude (and hence absolute i-b and magni- 
tude) is found, that does not contain a DLA. As I York et al.l 
l|2006l ) point out, this matched non-absorber QSO spectrum 
will be of very similar data quality and hence sensitivity 
to similar levels of absorption at each wavelength, because 
of the regulation of exposure levels directly at the telescope 
and the exposure level check of all SDSS spectra before they 
are admitted to the database. 

Given that there is a dispersion in the spectral slopes of 
QSO of similar magnitude and redshift, and hence a dif- 
ficulty in matching spectra one by one, the strength of the 
method proposed here relies upon averaging over an as large 
as possible set of QSOs with individually different slopes to 
obtain a robust measure of the (geometric) mean spectra for 
both the samples of QSOs containing DLAs and those free 
of DLAs. SDSS provides for both the largest homogeneous 
samples of DLAs as well as the largest reservoir of QSOs 
from which to draw comparison spectra. 

3.1 Constructing the non-absorber sample 

To find an appropriate match for a DLA-containing 
QSO spectrum, we search the SDSS DR5 databasqj 
after cleaning it from BALs and also the DLAs found 

^ Currently, the cataloque of I Gibson et al. I 1I2OO9I ') only allows 
us to derive BAL-free samples up to DR5. Hence, we have to 
restrict ourselves to finding matches in DR5 only, as a reliable 
BAL search is beyond the scope of this paper. We have, in fact, 
checked for the effects of possible BAL inclusion by comparing 
samples drawn from DR5 strictly to those drawn including 'pure' 
DR7 QSOs, and concluded that indeed they lead to substantially 
'redder' matching composites. Note tha t both the DR5 and DR7 
DLA samples rely upon a BAL removal llProchaska fc Wolfel2009l: 



in both lists. This resulted in a pool of 8371 QSOs 
with redshifts 2.24<Zem <5.31, from which to draw the 
comparison spectra. We minimise the difference of the 
match to the QSO in the two-dimensional parameter space 
of the i magnitud^l and emission redshift, defined as 
jj2 _ ( ^^q'^q7°"' Y' + ((H^)^- Note that we have corrected 
the i magnitudes in each case for Galactic reddening. The 
normalisation factors in the distance metric were chosen 
such that the number of outliers in both quantities are 
minimal, as detailed in the following paragraph. 

Figure [2] shows the distribution of deviations between 
the DLA containing spectra and their best matches in 
redshift and brightness for the two samples. The shape of 
the brightness distribution of the quasars in SDSS within 
the redshift range probed here, introduces a systematic 
Malmquist bias : in general there is a slightly higher chance 
to select a fainter match at higher redshift. This effect 
is particularly pronounced in the brightest DLA QSOs, 
because the number of potential counterparts below i mag 
^18.0 becomes too low to select an adequate match. Hence, 
we decided to exclude all DLA sightlines from the following 
analyses that deviate more than 0.08 in i mag and/or 0.007 
in relative redshift difference (dotted lines in the two panels 
of O. This reduces the DR5 sample to 475, and the DR7 
sample to 676 DLAs. The mean 1 band magnitude difference 
Ai(DLA - match) for these is -0.0029, and -0.0022 mag, 
whereas the mean relative redshift differences are A2:(QS0 
- match) /z = 9.3 x 10"^ and 5.3 x 10"^ 

In principle, we could increase the number of matches 
for each DLA spectrum to more than one counterpart. 
However, the gain in S/N (or reduced scatter) for the ratio 
of DL A/match, dominated then by the DLA composite's 
noise, which cannot be reduced, does not compensate for 
the loss in matching precision : the more matches are being 
culled from the pool, the more severe the bias becomes. 
Of special importance here are they increasing redshift 
differences that may distort the precision at which the 
composite match spectrum reproduces the shape of the 
Hence, we have decided to utilise only the best match in 
each case. 



3.2 Constructing composite spectra 

The redshift distribution for the DLA systems in conjunc- 
tion with the SDSS spectral coverage (3800 A 5C A 
95OOA) allows in principle the construction of a compos- 
ite spectrum in the DLAs' restframe from ~ 600 A to ~ 
3000 A. Obviously, however, the flux distribution of each 
QSO blueward of 1215 A in the quasar restframe is strongly 
suppressed due to the neutral H column densities sampled 
inside the forest, in a manner specific to each sightline dif- 
ferently. Hence, in order to reliably determine the effects of 

iNoterdaeme et al.ll2009| ) , and hence it is mandatory to implement 
such a check for BALs into the matching strategy, as well. 
^ We have used the values for the psf 'target' magnitudes and 
Galactic dust reddening as well as the QSO redshift SpecZ listed 
in the QSOConcordanceAU DR5 and DR7 datatables. 
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QSO i magnitude [mag] 

Figure 1. The distributions of the two DLA sample sighthnes in QSO emission redshift and the i band magnitudes. Plotted are also the 
distributions for the underlying SDSS DR5 and DR7 QSOs. Note the factor 2 increase of sources from DR5 to DR7. The decrease of DLA 
containing sightlines towards lower redshifts, despite an increasing number of QSOs, is due to the ever shorter pathlengths over which 
HI absorption can be identified for the spectroscopic coverage of the instrument trailing off at 3800 A on the blue side. In these, and the 
following figures, we have corrected the observed magnitudes for Galactic reddening. The insert shows the histogram of the absorbers in 
both samples. The majority of DLAs (80%) are located between 2.3 ^ z ^ 3.3. The shape of this distribution directly affects the length 
of the leverage arm in the DLA restframes that is at our disposal in the composite spectra. 



dust reddening by comparison with an unreddened contin- 
uum, we have to restrict ourselves in each spectrum to the 
part red ward of the beginning of the Lyman forest. We have 
therefore chosen to exclude all pixels in a spectrum that lie 
below 1240 A in the QSO restframe. This introduces a cutoff 
in the restframe of the absorbers that varies from spectrum 
to spectrum individually depending on the combination of 

Zem and Zabs- 

Geometric mean spectra were then generated from the ab- 
sorber and non-absorber sample in the following manner : 
we extend the SDSS wavelength grid with a stepsize of 0.001 
in log A to the lowest wavelengths covered by the spectra in 
the DLA restframe, and, as specified above, redward of the 
Lyman forest of the QSO itself For each wavelength data 
point Xdla.rf in this grid we compute for each sightline 
separately which pixel in the original spectrum comes to lie 
closest when transforming by division through the absorber 
redshift, i.e. we take the nearest integer value of 

N = ilog{\DLA,RFil.O + ZdLa) - log\start)/OmQl (1) 

where Xatart denotes the wavelength of the first pixel in 
each spectrum, zdla is the absorption redshift of the DLA, 
and A'^ is the pixel number starting from the blue end of the 



spectrum. Before computing the geometric mean of the fiux 
distribution in the sample of the thus determined pixels, 
we clear it of those individual pixels fl agged by the spe ctro- 
scopic pipeline as possibly corrupted (|Stoughtonll2002l ) . We 
have made no attempt in masking additional pixels, e.g. at 
the location of known absorption at different redshifts than 
the DLA. There is also no weighing of pixels according to 
e.g. S/N, as this possibly distorts the power-law nature of 
the combined spectrum. We ha ve corrected the indi vidual 
spectra for Galactic reddening (jSchlegel et al.l[l998l ). The 
geometric mean composite spectra for the non-absorber 
sets were generated in exactly the same way, using the 
redshift of the corresponding DLA in each case for equation 
[1] Note that we use in each case the spectra delivered 
by the DR7 pipeline, which includes the latest version of 
the spectrophotometric calibration. The advantage of the 
geometric mean is that it yields a meaningful average of 
the extinction at each wavelength by essentially averaging 
optical depths. 

Note that our method does not rely upon any inter- 
polation of the original spectra - at the price of sacrificing 
spectral resolution, because we bin together fiux values 
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Figure 2. The difference in i band magnitude and relative Zem emission redsliifts of the QSOs with DLAs and their matches for the 
DR5 (cyan) and DR7 (blue) samples. Note that there is a bias to pick fainter sources as matches, which is especially pronounced towards 
the brighter objects (i mag STIS.O), where the number of QSOs rapidly dwindles (cf. figure [TJ- The dotted lines indicate the location of 
the boundaries that we deem acceptable for a good match. 



that may reside in the DLA restframe separated by 1/2 
of a pixel at the extreme. We have tested the effects of 
this by comparison of the resuhing composite spectra to 
two cases of a more refined treatment : I. rebinning each 
spectrum after deredshifting it, and II. computing the 
weighted average of the fiux values of the two pixels in 
the original spectra that always straddle one wavelength 
element in the new DLA restframe grid. As expected, the 
differences over large wavelength ranges are negligible, and 
manifest themselves to a noticable degree only on a two 
to three pixel scale. Hence, we believe our method, being 
as close as possible to the data, is robust enough not to 
have to rely on the intricacies of interpolation as in case I. 
Because we have to throw out pixels with fiuxes f^ ^ 0.0, 
method II from above results in having to flag out twice as 
many datapoints, which is therefore deemed less acceptable 
despite the higher spectral resolution. 

Figure [3] shows the number of spectra contributing to 
each wavelength bin for the construction of the composite 
spectra in the DLA restframe. It also highlights to which 
location in these spectra the observed central wavelength 
of the i band (7481 A) is shifted in the DLA frames. 
Furthermore we have indicated the average brightness of 
the QSOs contributing at each wavelength, and the average 
DLA absorber redshift. Note that the region covered by 



the bulk of the spectra extends from ~1400 to ^--^2300 A, 
and exhibits fairly constant values of the mean brightness 
and absorber redshifts. The need for staying redward of 
the Lyman a emission, and the decreasing high (restframe) 
wavelength with increasing redshift of each QSO, result 
of course in fewer spectra contributing to the composite 
towards the edges of the wavelength regime covered by our 
samples. Hence, these spectral areas carry less weight due 
to their increased noise/scatter in the reddening analyses 
(as detailed in 3.3.). 



3.3 Derivation of extinction 

The geometric mean composite spectra of the 475 (DR5) 
and 676 (DR7) absorber spectra, and the corresponding 
non-absorbers are shown in the top panel of Figures |4] and 
[S] Note how the QSOs' emission lines for both samples are 
equally smeared out due to the different shifts applied for 
each specific combination of Ze,n and Zats , whereas only for 
the DLA sample the narrow absorption lines appear in the 
expected places. Because of the lower numbers of spectra 
contributing to each wavelength bin towards the blue and 
red end of the composite, the noise (or scatter) starts to 
increase. The rise of the mean flux towards the red side 
is an effect of the QSOs contributing to this part of the 



6 S. Frank and C. Peroux 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




1500 2000 2500 3000 

Absorber restframe wav elength [A] 

Figure 3. The number of spectra contributing to each wavelength bin for the construction of the composite spectra in the DLA restframe 
(sohd, thick hnes). Also shown is the histogram of the position of the central wavelength of the i band (7481 A), that is used to find 
matches (frequency scale units for the histogram on the right). The dotted line marks the average brightness of the QSO sample going 
into the composite at each position (see scale on right axis), and hence explains the rise towards long wavelength. This is accompanied 
by a similar drop in the average absorber redshift (dashed lines). 



spectra being on average of lower redshift and brighter 
(cf. Figure [ij . Additionally, the abrupt change in the flux 
beyond 2800 A can be attributed to the stack becoming es- 
sentially the average of a very few spectra from this point on. 

The lower panels of Figures |4] and [5] show the ratios of 
the two composite spectra. The cyan dotted line indicates 
the ratio expected assu ming an e xtinction of <E(B-V)> = 
0.01 for an SMC curve (|Peilll993 ), and represents th e uppe r 
limit on <E(B-V)> derived by iMurphv fc Liskd (|2004 ). 
It is immediately obvious that very little room is left for 
the signature of any dust extinction in the DLA sample. 
We can quantitatively assess how much extinction can be 
tolerated by our data with the following procedure. 

Suppose an unknown DLA stacked spectrum, 
/Di,A,nr(A) is reddened by dust within the absorbers, 
fDLA,r{X). Then, the ratio of these two spectra would sim- 
ply yield the extinction : ro(A) = foLA.nrW/ fDLA,r{>^) = 
IqOAAx ^ -pj^g underlying assumption for our method is that 
the stacked spectrum of the matches is a good representa- 
tion of the unreddened DLA stack. Our matching scheme, 
however, leads us to an unreddened spectrum fmatch,nrW, 
that, by definition, has the same flux as the (supposedly) 
dust reddened DLA spectrum at some fiducial wavelength 



Aq. Hence, the ratio of the two composites in this case is 
modified by the extinction ^4^0 at this matching wavelength 

tch.nr (A)//l3LA,r(A) 
= ;^q0.4x(A^-A^j,) ^2) 

Given an extinction law, the ratio becomes a function 
of E(B-V) only : ri(A) = ri{E{B - V). H ere, we have 
chosen the SMC extinction curve oflPeil l|l992h . and assume 
an Rv ~ 2.93 as typical for the total-to-selective ratio. Note 
that the choice of the SMC reddening curve is the most 
constraining regarding the absence of dust signatures - LMC 
or even MW dust type induces stronger reddening. Because 
we are fitting for the slope of this function, the potential 
offsets in Axg are of no relevance in this circumstance, and 
hence the normalisation of the ratio curve is not important. 

Figure [6] shows the results of these one-parameter 
fits for DL As (this work) and two SPSS samples of Mgll 
absorbers ijYork et al ] I2OO5I : IVanden Berk et al.l I2OO8I ). To 
obtain the latter two, we have applied exactly the same 
averaging scheme as outlined above. Before the calculation 
of the goodness of fit statistic, x^/'^ (where 1/ denotes 
the degrees of freedom for the fit), we have taken out 
areas of ±5 A from the known, narrow metal lines in the 
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Figure 4. Upper panel: The composite restframe spectra of the DLA DR5 sample (blue) and the corresponding sample of best matches 
(red). Before averaging the spectra, as detailed in the text, a Galactic dust extinction correction is applied. The stretched out broad 
emission lines have very similar shapes, indicating the suitability of the match for the analysis at hand. Only the DLA composite exhibits 
narrow metal line absorption features. The rise of the curves towards the red end is purely an effect of the larger average brightness of 
the sources at lower redshift, which are the sole contributors at this wavelength regime (cf. Figure[TJ. Lower panel: Ratio of the the two 
comp osite spectra from the upper panel. The cyan dotted line indicates the ratio expected for an SMC extinction with E(B-V)=0.01 ijPeil 
Il992f) . which exhibits clearly much more curvature than the data allow. The black dotted vertical line at 1700 A indicates the position 
bluewards of which the strong emission lines of CIV, possibly with contributions from Fe lines, and SilV render the exact matching of 
spectral shape increasingly difficult (see text for details). 



spectra. We find that we can constrain tlie amounts of 
extinction for tliose two samples of Mgll absorbers very 
tightly : <E(B-V)> = 0.012±0.002 for the intervening, 
and <E(B-V)> = 0.025±0.003 for the associated cases . 
We there fore reproduce the r esults of lYork et al.l (|2005l ) 
and IVand cn Berk ct alj (|2008l ) well, despite the fact that 
we do not explicitly also remove regions of the spectra 
that are affected by strong absorption at redshifts dif- 
ferent from the Mgll absorbers. We take this test as a 
sign that both our stacking algorithms and the method 
to estimate the reddening are robust and accurate for 
these Mgll samples. We note, however, that the lower 
redshift of these absorbers (compared to our samples of 
DLAs) may reduce the problem of the foreground red- 
dening due to the lower redshift path sampled, and the 
intrinsically lower density of absorbers in the low-z universe. 



We have also applied this spectral method to the 
samples of DLAs and 'pool' QSOs used by Vladilo et 
al. (2008) for their photometric analysis. These authors 
have taken care to take out QSOs that may be afflicted 



by absorption unrelated to the DLAs. Again, we are 
able to reproduce their results very well, and derive an 
< E{B - V) >= 6.2 ± 1.8 X 10"^, indicating that both 
methods yield similar results when applied to the same 
samples. Here, the strength of our method is specifically 
that all the unrelated effects will affect both the DLA and 
the pool spectra in the same way, and hence we choose 
to increase the sample size to ensure that these effects 
are averaged out. We are making use of both the largest 
databases for DLAs and comparison spectra (increasing 
the number of DLAs by almost a factor of 3 for the DR7 
sample, and the number of the objects in the comparison 
pool by a factor of 5 compared to the Vladilo samples). 



A word of caution is necessary when trying to quantify 
the confidence intervals for our results. Nominally, the 
number of degrees of freedom (d.o.f.) for the two samples 
are 2543 (DR7) and 2583 (DR7), after taking out the pixels 
affected by metal-line absorption, and restricting ourselves 
to the wavelength region below 2500 A to avoid the large 
drop in both numbers of spectra contributing and average 
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Figure 5. Same as figure |4] for the DR7 sample. 
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QSO brightness. However, as we have detailed in the sec- 
tion on the construction of the composite spectra, adjacent 
pixels are not truly independent of each other. We have 
tested the robustness of the fits by binning the composites 
(with bin sizes up to 5 pixels, hence reducing the d.o.f. by 
the same factor). Hence, we arrive at measurements of the 
mean <E(B-V)> =-0.0013± 0.0025 (DR5) and <E(B-V)> 
=-0.0017±0.0022 (DR7) for the full samplefl, where the 
error budget adopts a generous 3(7 level for the approach 
with the 5-pixel bin. 

In order to test whether we could have detected an 
<E(B -V)> distribution as predicted bvlPpntzen fc Pettinil 
l|2009h and measured by IVladilo et al.l |2008t ). we have 
artificially reddened each DLA quasar spectrum by a value 
drawn from a lognormal distribution of width 0.4 dex and 
mean of log < E{B — V) >= —2.2, broadly consisted with 
those two studies. Figure [7] shows that we not only could 
clearly detect such a reddening, but also could retrieve the 
input value very accurately (red curve in Figure |6]|. In fact, 



* Of course, negative reddening values are per se unphysical, but 
the natural outcome of the fitting routine. These are largely driven 
by the fact that the regions in the spectrum with the strong QSO 
emission lines between 1300 and 1700 A cannot be expected to be 
perfectly matching, given the complicated line emission structures 
of quasars. This difiiculty, and not sample size, probably also sets 
the limiting precision of our method. 



the two X curves for our spectral analyses and the fake 
reddening are excluded by several a. 

We have also performed a 'null test' by comparing the 
composite spectra of the matched samples with composites 
created from the samples of the 2nd best, 3rd best, etc.... 
matching samples. As expected, the analysis of the ratios of 
these stacked spectra yields no detectable reddening, and 
the range of <E(B-V)> values obtained by this procedure 
ranges from -0.002 ^<E(B-V)> s;0.002 for DR5, and 
-0.0018 ^<E(B-V)> ^0.001 for DR7. Hence we encompass 
the two measurements with the DLA samples. We note in 
passing that the main difficulty in matching the shape of 
the composites and their comparison counterparts lies in 
the strong emission lines of the QSOs, especially blueward 
of ^1700 A, where the Fe-complexes have a profound 
influence on the flux above the continuum. Indeed, when 
inspecting the ratio curves for the two samples, around 
that wavelength there seems to appear a 'kink' in them, 
which may be caused by a mismatch in the emission lines. 
This is the main driver towards these (very low) 'negative' 
<E(B-V)> values we obtain formally from the fitting. If we 
restrict ourselves to either the ratios blueward or redward 
of this feature, the fitting results in exact null values (at 
the expense of leverage in wavelength). 

One possible complication in matching the quasars' lu- 
minosities based upon their emission redshifts and the ob- 
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served i band magnitudes is that the latter could systemat- 
ically underestimate the true continuum flux for the DLA 
sample compared to their non-DLA counterparts because of 
a higher fraction of strong metal absorption lines due to the 
presence of the DLA. A rough estimate of the size of this 
effect indicates that in sightlines where a very strong ab- 
sorber (of the order of 5 A restframe equivalent width) is 
present in the i band, it will result in a photometric mea- 
surement that is decreased from the 'normal' QSO bright- 
ness by 2%. Thus we may in such cases be selecting QSOs 
as matches that are systematically too faint. We have tested 
whether such a possible bias may influence our analysis of 
the dust reddening signal by creating samples of matches 
following exactly the procedures described above, but in- 
troducing an artificial offset 5i in the i band magnitude by 
minimising A{iDLA + Si) — imatch- Even shifts as large as 
5i=-0.2 magnitudes do not change the overall shape of the 
composite spectra, yielding almost exactly the same results 
for the reddening estimates when fitting for the slope of the 
ratio curves. 



4 EXTINCTION ANALYSES FOR 
SUBSAMPLES 

While the full samples of DLA sightlines do not show any 
appreciable sign of dust reddening to a high degree of 
accuracy, it is still possible that there are certain subsets 
of DLAs and/or sightlines that do exhibit an amount 
of wavelength dependent extinction, which is simply too 
small to be detected when averaging over such large 
datasets. In order to check for such potential trends, 
we have divided the DLA lists into subsamples accord- 
ing to various criteria concerning the DLAs themselves 
(namely absorption redshift, HI column density, and the 
ab- or presence of metal lines) or the underlying QSO 
sightlines (apparent brightness or emission redshift). To 
keep a high enough number of DLAs we have in almost all 
cases decided to simply split the full samples into two halves. 

Figure [5] shows the results of a split (for the DR5 
samplejj into the highest (left panel) and lowest (right 
panel) HI column densities for the DLA. The mean of these 
subsamples are at log Nfr/=20.91, and log N/f/=20.42. 
Clearly, there is no difference in the fact that there is a null 
detection for both subsamples. 

Figure [TT] shows how differentiating the samples ac- 
cording to the absorber redshifts affects the dust reddening 
analysis. Again, both for the low and the high redshift 
subsamples (mean Za6s=2.63, and 3.39, respectively) there 
is no reddening detectable. Note how the higher redshift 
sample leads to a lower wavelength cut-off, and due to 
the generally fainter spectra a slightly higher noise in the 
composite. 

Figure 1101 demonstra tes that the DLAs classified by 
iProchaska fcWolfel l|2009l ) as containing no visible narrow 

^ We have performed the subsample analyses for both DR5 and 
DR7, yielding the same results, but only show the results for DR5 
here as examples. 



metal absorption lines in individual SDSS spectra, do also 
not exhibit any such signs in their composite. Only 88 of 
the 475 DLAs fall into this class (18.5%), and hence the 
noise in the composite is much higher than in the other 
subsamples. Yet, both classes do not show reddening. 

The separation into QSOs of low and high redshifts 
(shown in Figure [121 mean Za6s=2.31 and 3.78) is similar to 
the cut in QSO brightness (mean i band psf magnitude of 
i=18.62, and 19.48 mag, corrected for Galactic reddening), 
because of the trend for fainter sources with higher redshift. 
The subsample with the brightest QSOs (for DR5jf| is the 
only one that does show a positive value of the colour 
excess, but not signifcantly different from a null detection : 
<E(B-V)> =0.0012±0.0025. 

To summarise, none of these simple subselections do 
give rise to the assumption that there are significant trends 
in reddening regarding either the observational constraints 
(brightness and emission redshift of sources) or DLA prop- 
erties (absorption redshift, additional narrow metal line ab- 
sorption, and HI column density). 



5 SUMMARY AND CONCLUSIONS 

While our analysis technique reproduces well past re- 
sults f or other classes o f abs o rbers at d i fferen t red- 
shifts l|Vanden Berk et al.l [2008; 'York ct al.' '2006'), and 
DLAs using other m ethods (Murphy & Liskc 2004) and 
IVladilo et al.l |2008| ) where some extinction could be 
detected, we could no detect any significant extinction in 
samples of DLAs drawn from the Data Releases 5 and 
7 of the SDSS. These findings are based on the largest 
DLA sample ever studied, but also on the largest sample 
of quasars from which to draw a reference spectrum. The 
implication is that most of the DLAs observed in SDSS 
have a very small dust content. 

lYork et al.l l|2006l ) have analysed SDSS spectra with 
intervening strong Mgll absorbers, and detect clear dust 
reddening. More specifically, their subsamples with the 
highest Mgll equivalent width (subsamples 8 and 25 in 
their notation. Table A4) show the highest degrees of 
reddening (<E(B-V)> =0.032 and 0.085). These strong 
Mgll absorbers are then classified, based on the column 
densities of the accompanying metal lines, as akin to low 
redshift DLAs. We note however, that a direct measurement 
of the HI absorption in these absorbers was not possible 
due to the low wavelength cutoff of SDSS, and the mean 
HI column densities for the subsamples are derived by 
using an SMC relationship between <E(B-V)> and HI. 
I Rao et aXi (20i)g) also caution that not all high EW Mgl l 
absorbers are indeed DLAs, although lYork et al.1 (|2006D 
po int out tha t the s pecific requirement of the QSOs probed 
bv lRao et all (120061 ) with HST being UV-bright introduces 
a bias against red QSOs, that may have lead them to an 
underestimate of the DLA fraction. Hence, it is a priori 

^ The same cut in the DR7 sample does not show such a detec- 
tion. 
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Figure 6. Results of the onc-paramctcr fit for the dust reddening in various absorber types. Plotted are the P^r degree of freedom 
distributions obtained by fitting the ratio of the composite spectra with a one-parameter extinction curve (for details see tex t). The two 
curves in magenta and cyan sh ow the results for the samples of intervening and associated Mgll absorbers in SDSS spectra jYork et al.l 
I2OO5I : I Vanden Berk et alfe oOS'). We recover the extinction <E(B-V)> measured by these authors. The black and blue curves are a result 
of the fit to the DLA sample. Note that we can clearly rule out any signs of extinction <E{B-V)> > 0.002 with very high confidence. The 
curve in red shows the result for the DR5 sample that has been artificially reddened by a lognormal distribution with mean <E(B-V)> 
=6.3xl0~3, which is well recovered. 



not clear whether our results contradict these findings. 
Also note that ther e may be a redshift evolution effect. In 
fact, iMenard et al.l (|2008l ) have noted that the rest-frame 
colour excess E(B-V) for their Mgll samples shows a strong 
decline with increasing redshift (their figure 10, lower right 
panel). Thus, it is possible that the difference in the mean 
redshifts of ou r samples ( z qps ~3 .1) compared to the lower 
Zabs ~1.4 of lYork et al.l ( 2006l l may partly explain the 
seemingly discrepant results. 

Yet another strategy involving spectral analysis is 
fitting the QSO continuum spectra by a power-law, and 
checking for deviations in the slope of these for spec- 
tra with and w it hout DLAs. Employing this technique, 
iMurphv fc Lis ke' (2004') can place an upper limit for a 
sample of 70 DLAs from SDSS DR2 of <E(B-V)> =0.01. 
Note that we can rule out such a high value for our samples 
with a very high degree of confidence. 

In addition, there have been claims of firm detections 
of dust reddening t owards sightlines affi icted by intervening 
DLAs. Specifically, IVladilo et al] lj2008l l report a mean E(r- 
z)=27(±9)xl0"^ for a subsample of the same DR5 DLA 



sample that we use. When applying our method to these 
subsamples of DLAs and the subsample of DR5 'pool' QSOs 
for the constructing of t he comparison spectr um, we recover 
their results. Note that Iviadilo et al.l (|2008l ) have removed 
all quasars from their samples that may have additional ab- 
sorption unrelated to the DLA, reducing the number of ob- 
jects in each sample by a large fraction. We conclude that 
the differences do not reside in the methods, but the sample 
selection. We note that the reddening distribution of unde- 
sired interlopers probably has a tail with a small number of 
objects exhibiting much larger than the typical reddening 
expected from DLA systems. It may hence be that for such 
low frequency events, that cases of e.g. BAL QSOs, despite 
efforts to filter them out, enter the DLA samples more often 
than the comparison pool, and hence could wash out the sig- 
nal. Increasing the sample size reduces the possible impact 
of such 'special', very red QSOs. 

Our result is also at odds with expecta tions from the 
evolution of high-redshift galaxies (but see iPontzen et al.l 
(j2008h . It should be emphasised however that there are 
two separate issues: i) what is the dust content of the 
quasar absorbers we know of today and ii) which fraction of 
the quasar absorbers are missed because of the extinction 
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Figure 7. The result of artificially reddening the DR5 DLA sample spectra with each QSO subjected to an <E( B-V)> drawn r andom ally 
from a lognormal distribu tion of mean log <E(B-V)> =-2.2 and width 0.4 dex, resembling the values found bv lVladilo et aL I l l2008f) and 
iPontzen &: Pettinil ||2009| '). The reddening is not only clearly detected by our method, but also the input value is returned accurately. 
Labelling and underlying comparison spectrum just like in Figure [4l 



they induce on their background quasar such that it is 
not selected in the firs t place. Indeed, models of galactic 
chemical evolution (e.g. IPrantzos fc Boissierl (|2000l ) predict 
the existence of regions of high metallicity and column 
density which are mostly not detected in current DLA 
surveys (|Boisse et al.lll998l ). One way to reconcile these two 
findings is to invoke a dust-bias in the observations. These 
apparent discrepancies could indeed be explained by the 
dust content of the absorbers themselves: they could absorb 
enough light from the quasar to make the latter too faint 
to be included in current quasar samples, which barely 
go beyond magni tude g~20.2. In fact, an extension of the 
I Fall fc Pea (|l993f ) calculations taking in to account recent 
observations addresses both these issues: IVladilo fc Perouxl 
propose that whilst the dust content of the DLAs 
in current samples is not high, the missing fraction is 
possibly quite important ranging from 30% to 50% at 
z = 1.8 — 3.0. One way to to test such a hypothesis is to 
look at fainter objects than what can be probed with SDSS. 
Using a Bayesian analysis of these issues in lig h t of d ifferent 
techniques and datasets IPontzen fc Pettinil (|2009l ) posit 
that probably less than 10% of DLAs have been missed 
in current optical samples, and also derive an <E(B-V)> 
posterior probability distribution peaking strongly around 
log <E(B-V)> ~ -2.0. 



On the other hand, if the absence of the detection of 
extinction in DLAs is not due to an observational bias, then 
selecting in absorption implies that we are picking a popu- 
lation of high-redshift galaxie s which have prob a bly u nder- 
gone little s t ar for mation, e.g. IWeatherlev et all l|2005ll and 
IWolfe et all (|2006l '). This is difficult to reconcile with the fact 
that these system are gas-rich by construction. 
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Figure 8. Comparison of the geometric mean spectra and their ratios for the subsamples with the highest (left panels) and lowest 
(right) column densities log J^hi- Note that for this, and the following plots (with the exception of the metal- no metal distinction), each 
of the subsamples simply contains half the number of sightlines as the original full sample. 
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Figure 10. Comparison of the geometric mean spectra and their ratios for the subsamples with the DLA that do not exhibit sign of 
metal-line absorption (left panels, 88 DLAs in total) and the ones that do (right panels). Note the conspicuous absence of any sign 
of metal-line absorption in the DLA stack on the left. 
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Figure 11. Comparison of the geometric mean spectra and their ratios for the subsamplcs with the lowest (left panels) and highest 
(right) DLA absorber redshifts. Notice the lower cut-off wavelength for in the DLA restframe for absorbers with higher redshifts. Despite 
the fact that hence the 'kink' beyond ~ 2500A in the stacked spectrum of the low redshift subsample gets even more pronounced than 
in the overall sample, the matched spectrum still follows this quite accurately. 
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Figure 12. Comparison of the geometric mean spectra and tiieir ratios for tiie subsamples with the lowest (left panels) and highest 
(right) QSO emission redshifts. Note the rapidly increasing noise for the high redshift subsample towards the red part of the composite. 
The separation into QSOs of low and high redshifts is similar to the cut in QSO brightness. The subsample with the brightest QSOs (for 
DR5) is the only one that does show a positive value of the colour excess, but not signifcantly different from a null detection : <E(B-V)> 
=0.0012±0.0025. 



